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'^ \ Abstract 

C ■ 

Q ■ The effect of a moving defect particle for the one-dimensional partially asymmetric 

, ^, . simple exclusion process on a ring is considered. The current of the ordinary particles, 

the speed of the defect particle and the density profile of the ordinary particles 

^^ are calculated exactly. The phase diagram for the correlation length is identified. 

^r^ ■ As a byproduct, the average and the variance of the particle density of the one- 

OO I dimensional partially asymmetric simple exclusion process with open boundaries are 

^ \ also computed. 

a^. 

5^ 1 Introduction 



Recently, the one- dimensional asymmetric simple exclusion process (ASEP) has attracted 
much attention in various fields of science including mathematics, physics and biology |]l|- 



fi 

^ ' There are several reasons for this. First of all, the ASEP has a lot of applications to realistic 

O ■ systems. For instance, the ASEP can be considered as a simplest model of traffic flow 0-0 . 

It is also used to study the dynamics of interface P]. There are other applications as well. 

Second, the ASEP shows rich non-equilibrium behaviors such as the shock wave P^TTI, 

rS \ boundary induced phase transition [0, and the unusual dynamical scaling |1^,|14|. Since 

c^ ' there is no established framework for far-from equilibrium systems, the ASEP has played an 

important role in non-equilibrium statistical physics. Third, some properties of the ASEP 

can be studied exactly by the Bethe ansatz or the so-called matrix product ansatz. Since 

the effects of fluctuations become strong in low dimensions, mean-field analysis sometimes 

fails to give correct answers. Exact solutions give us an insight how to understand the 

problems which can not be solved exactly. 
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The ASEP is a one- dimensional lattice gas model. Each particle tends to hop to the 
right with rate p and to the left with rate q. In addition, they interact with one another 
through hard-core exclusion interaction. In this article, we consider the stationary state 
of the ASEP on a ring with one defect particle. We sometimes refer to the particles as 
ordinary particles to distinguish them from the defect particle. We assume that the defect 
particle tends to hop to the right with rate a and that it can exchange its position with 
an ordinary particle on the left nearest neighboring site with rate (3. By rescaling time, we 
can set p = 1 without loss of generality. Then the above rules are represented as 

10 -^ 01 with rate 1 

01 -^ 10 with rate q . ^ 

20 -^ 02 with rate a ^ ' ' 

12 -^ 21 with rate /3. 

Here 0, 1 and 2 denote an empty site, an ordinary particle and the defect particle re- 
spectively. Several exact results have already been obtained for the q = case |I^, p!6[. 
This case is referred to as the "totally asymmetric" case in the following. The current of 
ordinary particles, the speed of the defect particle, the correlation length and the current 
fluctuation have been computed exactly for this special case. Much less is known for the 
"partially asymmetric" case where q is not necessarily zero although some exact results 



were obtained in ^Tj. In this article, the model is studied for the < g < 1 case. We com- 
pute the current of the ordinary particles, the speed of the defect particle and the density 
proflle of the ordinary particles seen from the defect particle. We exploit the connection of 
the present model to the partially ASEP with open boundaries, which was recently solved 
by using the theory of the g-orthogonal polynomials |]18|,|19[. We do not consider the q > 1 
case, in which the current becomes zero in the thermodynamic limit. We remark that, for 
the partially ASEP with open boundaries, the q = I case and the g > 1 case were studied 



in ^0\ and ^1\ respectively 



The paper is organized as follows. In the next section, the stationary state of the 
process is constructed by the so-called matrix product ansatz. The current of the ordinary 
particles and the speed of the defect particle are calculated in section ^. In section H, 
the average density proflle is computed. The phase diagram for the correlation length 
is identifled. The concluding remarks are given in the last section. In Appendix A, we 
compute the average and the variance of the particle density for the partially ASEP with 
open boundaries. This is a byproduct of the computations in the main text. 

2 Exact Stationary State in Matrix Product Form 

The stationary state of the process can be constructed by applying the so-called matrix 
product ansatz. This technique was flrst introduced in [Q to study directed animals. 



It was then applied to the ASEP in p3[. Since then this method has been successfully 



generalized and applied to many interesting phenomena. See for instance the references 



in [|T8|. In the formalism of the matrix product ansatz, the stationary state of the pro- 
cess ( |1.1|) is constructed as follows. Let us denote the number of particles as N and 
the lattice length as L + 1. To be specific we explain the construction mainly by using 
an example for the case L = 2, A^ = 1: the case where there are one ordinary particle, 
one defect particle and one empty site on the chain with three sites. For this exam- 
ple, the stationary state can be written as a linear combination of the six configurations 
012, 021, 102, 120, 201, 210 and is completely determined by specifying the corresponding six 
coefficients P(012), P(021), P(102), P(120), P(201), P(210). Each coefficient P(ri, t^, rg) is 
the probability that the system is in a configuration T1T2T2 with Tj = 0,1,2 {j = 1, 2, 3) for 
the stationary state. 

Step 1: To each configuration 0, 1, 2 at one site, we associate a matrix £", D and A respec- 
tively. The space on which these matrices act is not specified at this stage. 

^ ^ 

1 ^ D (2.1) 

2 ^ A 

Step 2: To each configuration tiT2T^ of the whole lattice we associate a matrix product. 
For instance the configuration 120 is associated with the matrix product DAE. 

Step 3: We construct a state, of which each coefficient is given by the trace of the matrix 
product. For the present example, the six coefficients p(")(012), p(")(021), p(")(102), 
p(")(120),p(")(201),P("n210) are given by Tr(EP)A), Tr(EAP)), Tr(P)EA), Tr(P)AE), 
Tx{AED),Ti{ADE) respectively. Here we denote the coefficients as P'^"^(ri, r2, ts) be- 
cause the corresponding state is unnormalized. The coefficients of the normalized state 
P(''"i5 ^2; T3) is obtained by simply dividing the unnormalized ones by the normalization 
constant, 

Zi=3,^=i = Tr(P)EA) + Tr(P)AE) + Tr(EP'A) 

+ TiiEAD) + Tr(AP)E) + Tr(AEP)). (2.2) 

In the above prescription, there is a problem about the existence of trace. We will see that 
we can take a special form of the matrix A which ensures the existence of the trace. 

Step 4'- We can show that the constructed state is the stationary state of the process if the 
matrices D, E, A satisfy the conditions, 

DE-qED = CiD + E), 
(3D A = CA, 
aAE = (A. (2.3) 

Here ( is an arbitrary number and is taken to he ( = 1 — q hereafter. In addition, we set 
a = a/(l — q), P = /5/(l — q). The demonstration of the fact that the above algebraic 



relations are sufficient conditions for the stationarity of the process proceeds in almost the 



same manner as for the ASEP [^] and therefore is omitted. As we will see at the end of 
the section, we can take a matrix A of the form A = \V){W\, where \V) and {W\ are some 
vectors. Then we notice that the above algebraic relations (|2.3|) are exactly the same as 



those for the partially ASEP with open boundaries [IS, 23|. 



Step 5: Obviously, the construction of the stationary state above works for any choice 
of L, A^. Explicitly, the coefficients P(ri, r2, . . . , r^) of the normalized stationary state is 
given by 



1 ^ 



(2.4) 



i=i 



for the configuration T1T2 . . -Tl which satisfies Ylj=i T?' ~ ^- Here the normalization con- 
stant Zl at is 



^L,N 



E {W\l[[r,D + il-r,)E]\V). 



(2.5) 



ri,T2,. . . ,TL i — 1 



The summation is over Tj = 0, 1 for j = 1,2,... ,L with the condition Ylj^i Tj = N. Hence 
we have constructed the stationary state in matrix product form. 



In this article, we will compute the current of the ordinary particles, the speed of the 
defect particle and the density profile of the ordinary particles seen from the defect particle. 
These quantities are written in terms of the matrices as follows. 
The current of the ordinary particles is given by 

J = Prob(rj = 1, Tj+i = 0) - gProb(rj = 0, Tj+i = 1) + f3PT oh {tj = 1, r^+i = 2) 

= (1 - g) iJ^^I^zl^K _ ^-^^^-i.A^-i j _ (2.6) 

The speed of the defect particle has a similar expression, 

v = a Prob(rj = 2, r^+i = 0) - /5 Prob(rj = 1, r^+i = 2) 

= (l-g) ^^-^'^~^^-^'^'\ (2.7) 

Zl,n 

The average density of the ordinary particles at site j seen from the defect particle reads 



\""jl L,N 



E 



Tl , T2 , . . . , Tj , ■ • ■ , TL 
Tl + ■ • • + ft + ■ ■ • + Ti = Af - 1 



{W\ WinD + (1 - n)E]D W [nD + (i - n)E]\v), 

k=j+i 



k=l 



f2.8) 



where fj indicates that it is omitted. This is an unnormahzed quantity. The normahzed 

density is given by {nj)L,N = {njy^^j^/ZL,N- 

Here we note that the process has a particle-hole symmetry which reduces our compu- 
tations greatly. The process is invariant when the particles and holes are interchanged, the 
direction of particle hoppings are reversed, and a and (3 are exchanged at the same time. 
In other words, the process is invariant under the change, 

particle ^^ hole 

a^ p (2.9) 

site number j ^^ site number L — j + 1. 

Due to this symmetry, it is sufficient to obtain the density for the right half of the system, 
when we look at the system by setting the origin at the position of the defect particle. The 
density for the left half of the system is obtained by using the above symmetry as 

(nj)L,Ar(a, p) = 1- {nL-j+i)L,L-N{P, a), (2.10) 

where the dependence of {nj)L on the parameters a and f3 are explicitly indicated. 

So far, we have been considering the situation with a fixed number of particles. In other 
words, our discussions have been for the canonical ensemble. We are mainly interested in 
the physical quantities in the canonical ensemble. They can be directly compared with 
the results of the computer simulations because they are usually done for fixed particle 
numbers. However, it turns out that computations are much easier in the grand canonical 
ensemble. Since the meaning of chemical potential is unclear for our model, we define the 
grand canonical ensemble simply as a superposition of canonical ensembles. Let us intro- 
duce the fugacity C,"^ which is associated with the ordinary particles. Then the coefficients 
Pgce(ti, T2, . . . ,Tl) of the normalized stationary state for the grand canonical ensemble is 
compactly represented as 

1 ^ 

PGCE(ri, r2, . . . , tl) = -^J^^{W\ Hl^.'^'^ + (^ " ^M\V)- (2.11) 

Here the normalization constant Zl{C,) is a summation of the coefficients of the unnormal- 
ized states. In matrix notation, this is simply written as 

Zl{0 = {WlC^lV), (2.12) 

with 

C = ^'^D + E. (2.13) 

We took the fugacity not as C, but as ^'^. This is only for the notational simplicity in the 
formulae in the subsequent sections. But this causes an ambiguity; ^ and —C, gives the 
same fugacity C,"^. In the following we sometimes give explanation only for Re^ > 0. But it 



is clear how the Re ^ < case should be considered. One remark is in order. If we set ^ = 1 
in ( p.ll| ), Pgce(ti, T2, • • • , tl) gives exactly the stationary state of the partially ASEP with 
open boundary condition in which there are particle input (resp. output) at the left (resp. 



right) end of the system with rate a (resp. (3) |]T8|,|23|. In other words, the partially ASEP 
with open boundaries can be regarded as a superposition of the present one defect particle 
model with all numbers of particles A^ = 0, 1, 2, . . . ,L with equal weight. This fact allows 
us to compute the average and the variance of the density for the partially ASEP with 
open boundaries. This is done in the Appendix A. 

Before closing the section, we give some notations, definitions and an example of the 
representation of algebraic relations ( |2.3| ). First we introduce g-shifted factorial. 



{a;q)n = (1 -a)(l - aq){l 
{a;q)o = 1. 



aq 



aq 



,n-l^ 



(2.14a) 
(2.14b) 



We also define 



j=0 



(2.15) 



for |g| < 1. Since products of g-shifted factorials appear so often, we use the notations. 



(ai, a2, ■ ■ ■ , ak] q)oo = (ai; q)oo{a2; q)oo ■ ■ ■ {ak] q)o 
(ai, 02, ■ ■ ■ , ak, q)n = {ai, q)n{a2; q)n--- {ak, q)n- 

An example of the representation of the algebra ( |2.3| ) is given by 



(2.16) 
(2.17) 



D 



1 


VI- 


q 















1 




Vi- 


q' 















1 




Vi- 


g3 



, E = D^, A=\V){W\, (2.18a) 



/ 1 \ 

b 



{W\=K,{a\=^\ 1, 



V(^'^)i' V(?;?)5 



5 ■ ■ ' 15 



\V)=^i\b), = K 



\/{q\q)i 
b^ 

\f{q]qh 



\ 



, (2.18b) 



/ 



where the superscript T indicates the transposition. Notice that the special form of the 
matrix A ensures the existence of the trace. We introduced 



1 1-g 

a = -l + - = -l + -, 

a a 

6^-1 + U-l + l^. 



(2.19) 
(2.20) 



From now on, we restrict our attention to the case where a > and 6 > 0, i.e., the case 
where a < 1 — g and (3 < 1 — q. This is mainly for the simphcity of the discussions. We do 
not expect there appears another phase transition if we consider a < or 6 < case. The 
constant k in (|2.18b|) is takes as k^ = (a6; q)oo so that (H^|^) = 1. An important fact about 



these matrices is that they are related to the so-called g-orthogonal polynomials p^-p6|. 



It is this connection that allows us to compute the physical quantities exactly in the 
subsequent sections. The relationship between the matrix product ansatz for the partially 
ASEP and the theory of the g-orthogonal polynomials was first clarified in [18| and has 



been exploited to study the partially ASEP with open boundaries in ||19|,pi 



3 The Current 

In this section, we compute the current of the ordinary particles ( |2.6| ) and the speed of 



the defect particle (|2.7|) in the thermodynamic limit. We need to know the asymptotic 



expression of Zl^n when L, N —^ oo with p = N/L fixed. 

3.1 Exact Formula for Zl{^) 

First we give the exact expression of Zl{S,) in the form of contour integral. Since the 
derivation is almost the same as that in [llH] , the proof is omitted. The result is 



(g,a6;g)oo [ {z\z-^-q)^[{l + iz){l + iz-^)]'^ 
^^^^ 47r^ Jc ' {ai-^z,ai-^z-\hizMz-\q)oo ' ^ ' ^ 

The contour C of the integral above is such that it includes all poles of the type a^~^q^ and 
b^q'^ with k = 0, 1, 2, ■ ■ ■ whilest it excludes all poles of the type a~^^q~'^ and h^^^^^q^^ 
with /c = 0, 1, 2, ■ ■ ■ . When a^^^ < 1 and 6,^ < 1 hold, the contour is simply a unit circle 
and Zl{C,) reduces to the single integral on the real line. But when a^^^ > 1 or 6^ > 1 
there appear other terms as well. In general, when 

<"' > ar'q >■■■> aC^^'^' > 1 > ar'g"*"'+' > ■ • • , 

b^>b^q>---> biq^''' > 1 > b^q-'''^' >■■■, (3.2) 

where n^°'' and n^^' are some non-negative integers, we have 

Zl{0 = Zf\i) + 4"^(0 + Z'^^ii). (3.3) 



Explicit expressions for the above three terms read 

4'^^(0 = (g,a6;g)oo5^«;]"^[(l + ag^■)(l + eVV^■)]^ (3.5) 

Zf (0 = {q, ah- g)oo 5^«;f [(1 + ehq'){l + rlg-^■)]^ (3.6) 

i=o 



with 



^^' (g,a6,a-i6e;g)oo(g,a&~i^-2g;g)i(i-e-2a2)g^'a3^&^' ■ 



(6) {h-'i-^-qU{h^e.Cih-q),{l-eh'q'^) 



Off course when a^~^ < 1 and/or 6,^ < 1, we should omit Zj^ and/or Zi^ in (^ 

The asymptotic behavior of Z]^ (^) can be evaluated by applying the steepest descent 
method, whereas the asymptotic behaviors of Zf^ (^) and Z\^ (^) are simply given by the 
j = terms in the summation. We find 



(^^Q~^g)o 

{a-'he-,q)c 
ir'b-';q) 
{ah-^^-^-q\ 



z^L^io - /5-?p.^jr [(1 + «)(i + ev^)]^ (3.10) 

^?^(o ^ KTj.'^t [(1 + ^'^)(i + ^"')]^ (3-11) 



as L — > 00 when Re^ > 0. 



3.2 Asymptotic Formula for Zl^n 

In this subsection we study the asymptotic behavior of Zl^j^. First we notice 

L 



z^(0 = 5^^L,^e'^. (3.12) 



Af=0 

The average density of particles is given by 



L^KX) J_j Oq 



Basically we can expect that each value of the fugacity ^^ corresponds to each value of 
the density. Then the density-fugacity relation is invertible and the equivalence of the 
canonical and the grand canonical ensemble holds. For our model, however, there appears 
a situation where the equivalence of the ensembles fails. For this case, we have to go back 
to ( |3.12D and use the inversion of it, 

^--iX/^P. (3.14) 



where Ci is a contour which encircles the origin anti-clockwise. Notice that ( |3.14| ) is 



always true for a finite L since Zl{^) is a polynomial of degree 2N. This formula also has 
an advantage that it allows us to obtain not only the exponent of the asymptotic behavior 
of Zi^iy but also the prefactor of it. 

In order to find the relationship between the density and the fugacity and derive the 
asymptotic expression of Zl^j^i, we need to know, for a give value of ^, which term in (|3.3| ) 
gives the main contribution to Zi{C,). For notational convenience we define ^o = p/(l ~ p)- 

When a + P > 1 {ab < I ) 
• Case < C, < a 

The main contribution to Zl{^) comes from Zj^. The density-fugacity relation ( p.l3|) 



gives ^ = a/o^o (= ia)- In order for this value to be in the region < ^ < a, the 
density p should satisfy p < 1 — a < ^. The asymptotic behavior of Zf^'j^ can be 
obtained by applying the steepest descent method. 

Case a < ^ < h^^ 

The main contribution to Zl{^) comes from Z\^ . The density-fugacity relation ( |3.13| ) 
gives ^ = ^0- In order for this value to be in the region a < S, < b^^, the density p 
should satisfy 1 — a < p < {3. The asymptotic behavior of Zj^ jy can be obtained by 
applying the steepest descent method. 

Case ^ > 6"^ 

The main contribution to Zl(^) comes from Zj^ . The density-fugacity relation ( p. 13 ) 



gives ^ = \/^o/b (= ^b)- In order for this value to be in the region ^ > b ^, the density 
p should satisfy 1 — a < (3 < p. The asyi 
by applying the steepest descent method. 



p should satisfy 1 — a < P < p. The asymptotic behavior of Z\^ jy can be obtained 



When a + P <l {ab > 1 ) 

• Case < ^ < \/a/b 

The main contribution to Zi(^) comes from Z^" . The density-fugacity relation ( p.l3|) 
gives ^ = ^a- In order for this value to be in the region < ^ < ^/a/b, the density p 
should satisfy p < P < 1 — a. 

9 



Case ^ > \fajh 

The main contribution to Zl{^) comes from Zj^ . The density- fugacity relation ( |3.13|) 
gives ^ = ^h- In order for this value to be in the region ^ > ya/h^ the density p 
should satisfy 13 <1 — a < p. 

Case ^ = \/ ajh 

So far we have not seen the values of the fugacity which corresponds to the region 
/3 < p < 1 — a. It turns out that this is related to the fact that the density-fugacity 
relation ( p.l3| ) breaks down at ,^ = \falh. Notice that Z}°jy(,^) and Z\j^{^ in ( |3.10D 



and ( p.ll ) have a pole at this point. Physically, this is related to the existence of a 



shock. It might seem difficult to know the asymptotic behavior of ^l^at for this case. 



But this is not the case. Let us employ (|3.14|) and take the contour as a circle C^_ 



with the radius R. As long as L is finite, we can change the radius R as we wish. 
To know the asymptotic behavior of ^L,Ar, we want to use (p.9|)-( pTTlD and apply the 



^W/ 



steepest descent method. It turns out that the saddle point associated with Z\ [^ 
(resp. Z\j{^') is at ^ = ^^ (resp. ^ = ^b) and the contribution from Z\^ '{^ is smaller 
than the others. Hence, to apply the steepest descent method, we have to take the 
radius of the contour as R\ = ^a(> yojb) for Zi'^^) and as R2 = ^bi< yojb) for 
Z\^ {^). At first let us take i? as i?2- We have 



.J_f d Ub-'r';q)oo [{l + eb){l + b-')]^ 
^ 2m Jc,^ e {ab-'^-'; q)oo ^^ ' ^ ' 

To apply the steepest descent method to the first term, we have to modify Cr^ to Cr^ . 
Since the integrand of the first term has a pole at \/a/b, there appears a contribution 
from the pole when the contour is modified. More explicitly, we see 

1 f dUa-'e;q)oo[il + a)il + ea-')]'^ 



1 f d^ {a'^e-, q)oo [(1 + a)(l + ^^a-i)]^ {a-'b-'; g)^(l + a)^(l + b-')^ 



(3.16) 



2m J c,^ i ia-^be;q)oo ^'^ (g;g)oo(a6 



Now we can apply the steepest descent method to the two integrals. But it turns out 
that the contribution from the pole gives the asymptotic behavior of Zl,n for this 
case. 

Combining the above results, we see that there are four regions in the a-(3 plane, in 
each of which the asymptotic behavior of Zl^n in the thermodynamic limit has a different 
form. We refer to these four phases as Phase A, B, C and D. The asymptotic expression 
of Zi^AT in each phase is given by 

10 



Phase A: a < 1 — p, P > p 

Z L N c::^ Z T j^ c::^ —7—;. ^ 1= ; r- ('J--L') 

' (&eo;g)oo v^a^-^(l-a)^p^+V2(i-p)i-^+i/2 ^ ' 

Phase B:a>l— p, /3<p 



(a^o"'; g)oo 72^/3^(1 - /3)^-^p^+V2(i - p)i-^+i/2 

Phase C:a>l — p, /3>p 
^ _ ^(0) Q;/3(a + /3+l)(g,a&g;g)oo 1 

"'''" "''^"(« + P-iW-pP(eo-^ag,eo&g;g)L2-LV^(i-pr-^- ^- ^ 

Phase D:«<1 — p, /3<p 

Z,^^ (^-^-^^(""^"^'^)- L ^ (=zi^). (3.20) 

(l-«)(l-/?)(g;g)oo «^-^(l - «r/3^(l - /5)^-^ ' ^'^^ 



These are simple generahzation of the resuhs in [|15| for the totally asymmetric case 



3.3 J and v in the Thermodynamic Limit 

Now it is straight forward to calculate the current of the ordinary particles J and the speed 
of the defect particle v in the thermodynamic limit. The results are summarized as follows: 

• Phase A: J = (1 — g)p(l — p) and v = {\ — q){a — p). 

• Phase B: J = (1 - g)p(l - p) and t; = (1 - g)(l - /3 - p). 

• Phase C: J = (1 - g)p(l - p) and i; = (1 - g)(l - 2p). 

• Phase D: J = (1 - g)[p(a - /3) + ;3(1 - a)] and t; = (1 - g)(a - /3). 

4 Density Profile 

Now we turn to consider the average density profile. As in the previous section, we first 
calculate the density in the grand canonical ensemble. In the grand canonical ensemble, 
the (unnormalized) average density at site j is defined by 

i-i L 

kn,t\i) = E (^1 ui^^^'D + (1 - ^^•)^]^'^ n [^^'^'^ + (1 - ^^)E]\^)- 

Tl,T2,...,fj,--- ,TL k = l k=j + l 

(4.1) 
11 



To translate the results in the grand canonical ensemble into those in the canonical ensem- 
ble, we need the formula, 



<"^>--i/. 



(»j>i."K ) 



ii^^m^- (4.2) 



Now we explain the main idea how to compute {nj)}^ {^). Details of the computation 
are not presented here since they are similar to those in |]19| . Since it is easier to calculate 
the density difference than the density itself, we rewrite the density at site j as 



L-l 



\^i 



)l = E((^^)i"^(0 - K+i)irHO) + K)i"^(0- (4.3) 

k=j 



At the right boundary, we have 

M^l\0 = j{W\C^-'\V). (4.4) 

Similar formula can be found for the left boundary as well. As for the computation of 
{nk)L\0 - {nk+i)L\0^ we notice that 

= ^^{W\C''^\DC - CD)C^-''-^\V) 

= ^^{W\C''-\DE - ED)C^-''-^\V). (4.5) 

By using the two facts that DE — ED is a simple diagonal matrix and that the resultant 
series can be summed up by using a formula for the g-Hermite polynomials, it turns out 
to be possible to derive an integral expression of {uk/^yO ~ ('^fc+i)! (0- If "^^ sum up 
these terms from k = j to k = L — 1, after some computation, we get 

L-l 

Y.i(''^)'^r\o - K+i)Sr^(0) = h+ h, (4.6) 

k=j 



where 



^r 1 \ r \3 f '^^1 f '^-^2 



h = 7(a6;g)oo(g;g)c 



J Zmzi J 27nz2 

{zi ^^^ zi zf; g)oo[(i + e^i)(i + e^r')]^-^ 



X 



« ^^1, < ^% \ qziZ2, qz^ ^^2 ^ ?^i^2 \ q^i ^^2 ^ ^^^2, h^z^ ^; q)oo{z2 + z^^ - Zi - z^ 



-In ' 



^/ L \ / n3 /" '^'^1 /" '^'^2 



/2 = -(a6;g)oo(g;g)c 



4 y Zt:%z\ j l'KiZ2 

1 ) ^1 ) ^2 5 ^2 



X 



^r,^, ^;^o,^o-g)oo[(l + e^l)(l + e^^')p■-M(l + e^2)(l + e^2'')]''"'■ 



(a^-i^i, a^-i^;^ \ gzi^a, g^i ^^2 ^ ^^1^2 ^ ?% ^^2 ^ ^^^2, ^^2 ^\ q)oo{z2 + z^^ - Zi - z-^^) 

(4.7) 
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When < a^~^, &^ < 1, the contours of Zi and Z2 are both unit circles. For other values of 
the parameters, the contours are modified so that the analyticity of Ji and I2 are ensured. 
Combining the results for (?t-j)^ (0 with the formula ( |4.2| ), we obtain the density profile 
of the ordinary particles in the canonical ensemble. The density at the right boundary is 
computed by using ( [4.4| ). The integral Ji gives us the density at bulk region, which turns 
out to be p except for the phase D. Off course this is consistent with the fact that we 
are now dealing with the system with the fixed density p. The integral I2 contains the 
information about how the density decays near the right boundary. The density near the 
left boundary can be found by noticing the symmetry (|2.1CI| ). The results are summarized 
in the following. 



Phase A (a < 1 — p and P > p; a > ^q and b < I/.C0) 

The average density in the bulk region takes the constant value p. The densities at bound- 
aries are given by 



{ni)L = P, {nL)i 



1 — a 



P 



-P- 



(4. 



Notice that the density near the left boundary does not change and takes the constant 
value p. Computation of the density near the right boundary shows that the phase A 
subdivides into three phases. 

• Phase Ai 

This phase corresponds to 



p<(3< 



P 



p+(l-p)g' 



P< 



\/(l -a)p 



^a{l- p) + ^{l-a)p 



(4.9) 



in the a-/3 plane and to 



q/Co <b< 1/^0, b > 1/v^ 



(4.10) 



in the a-b plane. Let us denote the distance from the right boundary as / = L — j + 1 
from now on. Then the density near the right boundary decays exponentially as 



(%)l = P 



{ab%-l){q,a-'b-'^,\,q) 



[I + afo^o) (1 + b) (a-ife-ig, b'^^^'q; q), 
with the correlation length. 



-exp[— //r] 



(4.11) 



-In 



«(1-P) ^ (1-")P' 



1-/3 



/3 



(4.12) 
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Phase A2 

This phase corresponds to 



<a<l — p, p> = ^ (4-13) 



p + (l-p)g2 ' ^a(l-p) + ^{l-a)p 

in the a-/? plane and to 

^0 < a < eog"', & < 1/v^ (4-14) 

in the a-h plane. The density profile near the right boundary decays exponentially as 



I ^ _ v^(l + <o)(Q^^^o;g)oo(g;g)L exp[-//r] 



with the correlation length, 



r 



-1 



-21n[v/«(l-p) + V(l-a)p]. (4.16) 



The decay of the density is not purely exponential but with algebraic corrections. 

• Phase A3 

This phase corresponds to 

a< ^'7/K. ^ /3> ' , (4.17) 

P+(l-p)g2' ^ p+(l_p)g 

in the a-/5 plane and to 

a>^oq-^ b<q^^' (4.18) 

in the a-b plane. The density profile near the right boundary decays exponentially as 

with the correlation length, 

r-^ = In- -4-, ; TT- (4.20) 

la + {l-a)q][p + {l-p)q] ^ > 

Phase B (a > 1 — p and P < p; a < ^q and b > l/^o) 

This phase is symmetric to the Phase A. All results can be obtained by using the symmetry 
(p.9|). The average density in the bulk region takes the constant value p and does not change 
near the right boundary. The densities at boundaries are given by 

(ni)i = l-i^(l-p), {nL)L = p. (4.21) 

a 

14 



Phase Bi 

This phase corresponds to 

1 — p < a < 
in the a-f3 plane and to 



1-p 



1- p + pq 



a < 



(i-m-p) 



(4.22) 



^og < a < ^0, b> — 



(3p+J{l-P){l-p) 



6 



(4.23) 



in the a-b plane. The density near the left boundary decays exponentially with the 
correlation length, 



-In 



Pp I (l-/^)(l-p) 



1 — q; 



a 



(4.24) 



Phase B2 

This phase corresponds to 



pq 



I- p + pq^ 
in the a-fi plane and to 



< a < 1 — p, a> 



{i-m-p) 



(4.25) 



i^'<h<i^'q-\ a<b<^ 



Pp+J{l-f3){l-p) 



eo 



(4.26) 



in the a-b plane. The density profile near the left boundary decays exponentially 
with the correlation length, 



r"^ = -21n 



(3p+J(l-P)(l-p) 



(4.27) 



Phase B3 

This phase corresponds to 



P< 



pq 



in the a-P plane and to 



1 - p + pq^ 



a > 



1-p 



1- p + pq 



b> ^Q q , a<^oq 



(4.28) 



(4.29) 



15 



in the a-b plane. The density profile near the left boundary decays exponentially 
with the correlation length, 

r"^ = In^ J^ . (4.30) 

[/?+(! -/?)g][l-p + pg] 

Phase C (1 — a < p < /?; a < ^0 and b < ^q^) 

The average density in the bulk region takes the constant value p. The densities at bound- 
aries are given by 

K)L = l-4(l-pr, {nL)L = i. (4.31) 

a (3 

The average density decays near the right boundary as 



(n,), = p-y^ii-^l. (4.32) 

The density decays algebraically and hence the correlation length is infinite. The density 
decay at the left boundary can be obtained from the symmetry relation ( |2.1CI| ) . This is the 



same as the result for the g = case in |15 | 



--i\ 



Phase D (a < 1 — p and f3 < p; a > ^q and b > ^^ j 

In this phase, there is a shock. Hence the average density at bulk region is not a constant. 

For the densities at the boundaries, we have 

(ni)i = P, {nL)L = l-a. (4.33) 

There are two regions in the system. One is the low density region with the density j3 at 
left side and the other is the high density region with the density 1 — a at right side. The 
two regions are separated by a sharp interface, which we call a shock. When L is large, the 
average density profile is described by a continuous function n{x) in terms of the rescaled 
variable x = j/L (0 < a; < 1). The derivative of n{x) is found to be 

^^~Cexp[L-F(x)], (4.34) 

where 

F{x) = xln[(l + a)(l + a-^G{x))] + (1 - x)ln[(l + bG{x)){l + b~^)], (4.35a) 

^ ^ -X + p - abjl - X - p) + ^y[{x - py + ab{l - x - p)]^ + 4a&p(l - p) 

26(1 - p) 

The constant C in ( [4.34| ) is determined by the normalization condition, 

"^ dn{x) 



Q (XX 



-dx = l-a-l3. (4.36) 
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It turns out that the function F{x) takes the maximum value at Xq = (1 — « — p)/(l— a— /5). 
If we expand F{x) around x = Xq to the second order and approximate ( [4.34|) by a Gaussian, 
we find 



dn(x) / L , ~ nxi I Lil — a — l3Y{x — Xq) 

-■ ' '1 — a — fj) exp 



dx \j 2'n[p{a - 13) + ^{1 - a)] \ 2[p{a - {3) + (3{l - a)] 

(4.37) 



This is essentially the same expression as that for the totally asymmetric case [T^. But 
it should be remarked that ([4.37|) is only a consequence of the above approximation; the 
correct formula is given by ([4.34|) with ( [4.35a| ) and ( [4.35b| ). 

We thus obtain the phase diagram shown in Fig. 1. As we have seen, the phase A (resp. 
B) for the current subdivides into three phases Ai, A2 and A3 (resp. Bi, B2 and B3). 



This phase diagram has a richer structure than that for the totally asymmetric case |]T5| 
for which the phase A (resp. B) for the current subdivides into only two phases and the 
phases A3 and B3 are not observed. The situation is analogous to that for the ASEP with 
open boundaries ]TB|,p 



5 Concluding Remarks 

In this article, we have considered the partially asymmetric simple exclusion process on a 
ring with one defect particle. We have computed the current of the ordinary particles, the 
speed of the defect particle and the density profile of the ordinary particles seen from the 
defect particle. The main result of this article is the identification of the phase diagram for 
the correlation length shown in Fig. 1. It turns out that the phase diagram has a richer 
structure than that for the totally asymmetric case. 

There are several possible applications and generalizations of the analysis in this article. 
First there are several models for which the partially asymmetric case has not been solved. 
For instance we can generalize the analysis in [§ for the totally asymmetric case to the 
partially asymmetric case. Second, we didn't study the case where g > 1 in this article 
since our main interest was the phase structure of the model and it is obvious that there 
is a phase transition at g = 1. But it would be interesting to apply the similar analysis to 
this case and compute the physical quantities exactly. In addition, if we consider the limit 
g ^ 1 carefully, the crossover behavior should be observed since this limit is associated with 
the change of the universality of the model from KPZ universality class to EW universality 
class 0]. Here we only remark that the corresponding ASEP with open boundaries have 



been considered in [^| for g = 1 case and in |21| for g > 1 case. Lastly it would be 



also interesting to generalize our analysis to the multi-species models ||27|-pl|. Compared 
to the ASEP, much less is known about these models. Several investigations are now in 



progress [32|. The results about these will be reported elsewhere. 
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Appendix A Density Fluctuation for the Partially 

ASEP with Open Boundaries 

In the main text of this article, we are mainly interested in the system with a fixed number 
of particles. The physical quantities are computed for a given value of the density p. When 
we consider the ASEP with open boundaries, however, the situation becomes different. The 
system is connected to particle reservoirs at boundaries; the total number of particles can 
fluctuate. As pointed out in the main text, the stationary state of the partially ASEP with 
open boundaries is obtained by setting ^ = 1 in ( 2.11|) . Moreover the probability that the 
partially ASEP with open boundaries has A^ particles is given by Zlj~^/Zl{^ = 1). By 



Af) 



using the asymptotic expressions of Zl(^ = 1) in [J18l and those of Z^'^'j^ with / = 0, a, 6, rf 
in ( p.l7D -( p^ ), it is possible to obtain the asymptotic behavior of the probability measure 
of the particle density when L — > cxd. In particular, we can compute the average and the 
variance of the density for the partially ASEP with open boundaries. In this appendix, we 
do not explain even the basic properties of the model and use some terminologies without 
definition. See [|T^,p!9|, ^|33| -p5[] for instance. 

First we study how the asymptotic expression of Z^^^ in ( p.l7|) - (|3.20|) looks like as a 



function of < p < 1. As an example we consider Zj^jy. Let us rewrite (|3.17|) as 



^L,Ar — 



(a ^^o;q)c 



y/2np{l-p){b^o;q)c 



-exp[L-S(p)], 



with 



Then it is easy to see 



S(p) = -In [a^-P{l - a)V(l - pY'"] ■ 



dS _ 
dp 



«(1-P) 
_{l-a)p 



(A.l) 



(A.2) 



(A.3) 



becomes zero at p = a. In addition, we have 

d^S 1 



(9p2 p(l _ p) 



;<o) 



(A.4) 
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for < p < 1. Hence ^'j^jv takes a single maximum value at A^ = aL. If we expand S{p) 
around p = a and substitute it into (|A.1| ), we find 

L 






a(l — a) 



exp 



^^^-^)" (A.5) 



2a(l -a) 



when p r^ a. The width of the peak if of order L ^/^. As L goes to infinity, we expect 
that the probability measure of the density is given by a Dirac measure Sa- Similarly Z\^ j^ 

takes a maximum value at p = 1 — P; Zj^j^ takes a maximum value at p = 1/2. Now we 
explain the asymptotic behavior of the probability measure of the density and compute 
the average and the variance of the density. 

• low density phase (a < |, /3 > a) 

From the results in the main text, we see the following. For small densities which 
satisfy p < 1 — a and p < /3, the asymptotic behavior of Zl,n is determined by Z^^j^. 

For larger values of p, we have Zl^n — Z)^ j^ when 1 — a < p < (3 and Z^^j^ ~ Z]^ \^ 
when p > 1 — a. As L goes to infinity, however, the main contribution comes from 
Z^j^ in this phase. 

By using ( |A.5|) , we compute the average and the variance of the density to get 

(p) = «, (A.6) 

((P-(P))^)-^^^. (A.7) 

The variance of the density goes to zero as L ^ cxd. 

• high density phase {f3 < ^, a > $) 

For large densities which satisfy p > 1 — a and p > (3, the asymptotic behavior of 
Zl^n is determined by Zi'j^. For smaller values of p, we have Zl,n — -Z^ltv "when 
1 — a < p < P and Z^^n — -^l at when p < 1 — a and p < /3. As L goes to infinity, 
however, the main contribution comes from Zj^j^^- in this phase. 
The average and the variance of the density are found to be 

(p) = l-/3, (A.8) 

((P-(P))^)-^^^. (A.9) 



maximal current phase (a > |, /3 > |) 

For small densities which satisfy p > 1 — a, we have Zl^n — Z^j^; For large densities 
which satisfy p > (3, we have Zl^n' — Zj^ j^. In the middle region of densities which 
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satisfy 1 — a < p < /3, the asymptotic behavior of Zl^n is determined by Z\^j^. As L 
goes to infinity, the main contribution comes from Zj^ j^ in this phase. 
The average and the variance of the density are found to be 

(p) = ^, (A.IO) 

{{P-{P)?)^^- (A.ll) 



This formula agrees with the one derived in |^ for the special case q = 0,a = (3 = 1. 
coexistence line (a = /3 < |) 

For small densities which satisfy p < a, we have Zl,n — Zi'j^; For large densities 
which satisfy p > 1 — a, we have Z^^j^j ~ Zj^ ^. In the middle region of densities which 
satisfy a < p < 1 — a, the asymptotic behavior of Zl^n is determined by Z)^ jy- When 



we set a = j3 in ( p.20[ ), we notice that Z)^ j^ does not depend on p. This implies that 



the probability measure of the density does not have a peak. In fact the probability 
measure w{p)dp in the thermodynamic limit is 

a < p < 1 — a 



w{p) = 1-2'^ ^ . (A.12) 

U otherwise 



Hence we obtain 



(P) = l^ (A.13) 

((p-(p)f)-m^. (A.14) 

Notice that the variance of the density is finite for this case. This is related to the 
fact that we see a shock on the coexistence line. The position of the shock can be 



anywhere |jTI],^,|3^. This leads to a macroscopic change of the total number of 



particles or the finite value of the variance of the density. 
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Figure Captions 



Fig. 1 : Phase diagram for the current and the correlation length. The phase boundaries 
for the current are represented by the thick solid lines; those for the correlation length are 
represented by the thin solid lines. The dotted lines are not the phase boundaries. They 
are drawn for convenience. 
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